QUANTUM TOMOGRAPHY, PHASE SPACE OBSERVABLES, AND 
GENERALIZED MARKOV KERNELS 



JUHA-PEKKA PELLONPAA 

Abstract. We construct a generalized Markov kernel which transforms the observable asso- 
ciated with the homodyne tomography into a covariant phase space observable with a regular 
kernel state. Illustrative examples are given in the cases of a 'Schrodinger cat' kernel state and 
the Cahill-Glauber s-parametrized distributions. Also we consider an example of a kernel state 
when the generalized Markov kernel cannot be constructed. 

PACS numbers: 03.65.-w, 03.67.-a, 42.50.-p 



1. Introduction 

Quantum homodyne tomography and 8-port homodyne detection are important tools in 
quantum optics for reconstruction of a quantum state of hght (see, |14j and references therein). 
In both measurements, there are associated observables, normalized positive operator measures 
(POMs), which give the measurement outcome statistics of the input state. In the 8-port 
(or double) detection, the corresponding observable is a covariant phase space observable 
generated by a kernel state K (see, e.g. (ID])- Recently, Albini et al. [2] found a POM Eht 
associated with the quantum homodyne tomography (see also [llj). 

The aim of this paper is to find a connection between the two POMs Eht and Ek when K is 
regular enough. We will show that is connected to Eht via a generalized Markov kernel. A 
drawback is that the generalized Markov kernel is not a positive function, so that, it cannot be 
interpreted as a transition probability from Eht to E^^-. However, E^ is a postprocessing of Eht 
and thus subordinate to the tomographic measurement Eht- 

This work is a direct continuation of our earlier work [13j where we considered only the 
Husimi Q-function case K = |0)(0|. We also showed that there is no hope to represent Eht as 
a postprocessing of E|o)(o| at least by using generalized Markov kernels. 

The structure of this article is the following: In section 2, we introduce the basic notations, 

definitions and some well-known results. We also recall the structures of the POMs Eht and E^ 

associated with the tomographic and 8-port homodyne measurements. We will use the Radon 
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and Hilbert transformations which are probably famihar to readers so that their properties 
needed in this article are collected in the Appendix. 

Section 3 contains the main results of this paper and some illustrative examples. First we 
note that the tomographic data is somehow 'sharper' than the 8-port data for the given input 
state T. Then we construct the generalized Markov kernels for POMs (with suitable kernel 
states K) and then expand them to Hermitc and Maclaurin series. Both series are used in 
the examples. The first example shows the structure of the kernel associated to the phase 
space observable generated by a 'Schrodinger cat' state. In the second example, kernels are 
constructed for Cahill-Glauber s-parametrized distributions (or more precisely, for operator 
measures associated to them). We will see that although the Husimi Q-distribution has a very 
simple kernel, the Wigner function has no kernel at all. Finally, we see that there does not 
necessarily exist a generalized Markov kernel for an arbitrary kernel state K. 

2. Notations and definitions 

For any measure space (fi,//), where // is a positive measure on Q, we let L^(Q) (resp. 
L^(Q)) denote the space of equivalence classes of integrable (resp. square integrable) functions 
Q ^ C. In the case of the Hilbert space L^(n), we let ( ■ | ■ ) I'^^n) denote its innerproduct which 
is assumed to be linear with respect to the second argument. When Q = M", n = 1,2, the 
measure /i is always the n-dimensional Lebesgue measure. We let S (M") be the Schwartz space 
of rapidly decreasing smooth complex functions on R". We identify »S(]R") with a subset of 
Li(R") nL2(R"). 

Let be a complex Hilbert space with an innerproduct ( • | • ) and an orthonormal (number) 
basis {|n) | n e N} where N := {0,1,2,...}. Let C{n) (resp. T{n)) be the set of bounded 
operators (resp. trace-class operators) on Ti. The set of states (density matrices), that is, the 
positive operators T G T(H) of trace 1, is denoted by T{H)^ . For any B G C{7i), we denote 
Bmn '■— {m\B\n), m, n so that B — Ylmn=o-^mn\'m'){n\ where the double series converges 
with respect to the weak operator topology. 

Define the lowering, raising, and number operators 

oo oo oo 

A Vn + + 1|, A* + l\n + l){n\, N A* A = 

n=0 n=0 n=0 

respectively. Physically the Hilbert space H and the above operators are associated to a single- 
mode optical field. We will, without explicit indication, use the coordinate representation, in 
which Ti, is represented as L^(R) via the unitary map |n) i— > where /i„ G 5(R) is the nth 
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Hermite function, 



and if„ is the nth Hermite polynomial]^ Recall that it is customary to denote (formally) 
%Ij{x) = {x\ip) and \x){x\dx = I (the identity operator), that is, {ipl^p) = J^{ip\x){x\ip)dx = 
J^tp{x)(p{x)dx for aWtp, ipen = L'^{R). 



Define^ the position operator Q := -^{A* + ^) and the momentum operator P := :^(^* — A), 
which, in the coordinate representation are the usual multiplication and differentiation opera- 
tors, respectively: {Qi/j){x) = xiIj{x) and {PiIj){x) = —i^{x). For any 6 G [0,27r), 

Qe := (cose)g + (sin^)P 

is the (self-adjoint) rotated quadrature operator and Qe : B(M.) — > C(H) is the spectral measure 
of Qe; we denote the Borel cx-algebra of any topological space fl by B{fl). It is easy to see that 
Qe{R) = e''^^Qo{R)e-'^^ for all R G B{R) and 6 eR. 

Define the Weyl operator (or the displacement operator of the complex plane )|^ D{z) = 
^zA*-zA which 

D{zy = D{zy^ = D{-z) and e''"" D{z)e-"''' = D{ze"') 
where z G C and ^ G M. In addition, for all Zi^ Z2 & C, 
(2.1) D{zi + Z2) = e-'^'^'^''^^D{zi)D{z2). 

Let q, p eR and z = {q + ip)/ V^- Then 

D{q,p) := D{z) = e'^^''"^ = e-^^P/^gipQ^-igP _ ^igp/2^-igP^ipQ ^ 

that is, for all G H = L^{R), (e*P«^/')(x) = e'P'^tpix), {e''^^^){x) = tp^x + q), and 

{D{q,p)i;){x) = e-^'P/^e'P^ijix - q). 



Using the polar coordinate parametrization z = pe*^, p = \z\ = \J \{q'^ + p"^) , = aigz 
arctan(p/g), that is, q = ^/2pcos9, p = ^/2psm9, one gets 

^2 2) ^iV2pQe — ^i(V2p cos e)Q+i{V2p sin e)P _ ^iqQ+ipP _ D(^—p 



^Hermite polynomials are given by the recursion relation Hq{x) = 1, Hi(x) = 2x, and Hn+i{x) — 2xHn{x) — 

2niJ„_i(a;), or by the Rodrigue's formula Hn{x) ~ ( — l)"e^^(i"e~^^/dx". 

means the closure of a linear operator B and z is a complex conjugate of z G C. 
■^Recall that Weyl operators are associated to a unitary representation of the Heisenberg group H, or to a 

projective representation of C = M^. 
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Moreover, the two-dimensional Lebesgue measure is 

(i^z = pdpdO = -dp^dO = -dqdp. 

2.1. Covariant phase space observables and Wigner functions. Let i^' be a trace-class 
operator on H. Define an operator measure Ek '■ B{C) 'C,(H) as follows: for all Z G B{C), 

Ek{Z):=- [ D{z)KD{zyd^z = ^ [ D{q,p)KD{q,prdqdp. 

Jz ^TT Jz 

Since Ei^(C) = {ti K)I we see that Ek is normalized if and only if tri^ = 1. Moreover, Ek 
is positive if and only if i^' is a positive operator. If K E T(7i)^ we say that the normalized 
positive operator measure (POM) Ek is a covariant phase space observable generated by the 
kernel K. In principle, any covariant phase space observable can be measured by using double 
homodyne detection (see, e.g. [TU] and references therein). 
The parity operator is 

oo 

n := R{7t) = ^(-l)"|n)(r2| = n+ - H" 

n=0 

where n"*" := J2T=o |2A;)(2A;| and n~ := YlT=o 1^^ + l)(2/c + 1| are projections onto the closed 
subspaces H'^ := II+Ti and H~ := II"?-^ (i.e., H = © 7Y"). Elements of (resp. H~) are 
even (resp. odd) functions. For any ip Eli. = L^(]R) we can write 

(n±v^)(x) = ^[^(x)±^(-x)] 

so that (Ilip){x) = ip{—x) or, formally, 11 = Jj^ \x){—x\dx. 

The Wigner function : C ^ C of a trace-class operator T is 



(2.3) 



W^{z') = W^{q',p') := -ti [TD{z')UD{zy] = -tr [TD{2z')U] 

TT TT 



TP 1 p27v poo 

— tT[TD{-p,q)]e-"^''^-'P'Pdqdp= — / tr [Te*^^'3eig~iV2p(^'cose+p'sme)^^2^^ 
47r^ J^2 Atc^ Jo Jq 



where z' = {q' + ip')/V2 E C (see, equation (2.2)). Thus, the probability measures R t— > 
tr [TQeiK)], 6 E [0, 27r), fully determine the Wigner function of any T E T{H)^. If is 
integrabl^ then 



/ W^{q,p)dqdp = 2 I {z)d^ z = ii [T]. 



'^This need not hold even for all pure states T — \rj){rj\ (see example 3). 
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Recall that the Wigner function is always square integrable and 

(2.4) ii[TK] = 2n I W^{q,p)W^{q,p)dqdp = 47i [ W^{z)W^{z)(fz. 

In addition, 

W^{q,p) = - [ f{q-x,q + x)e^'P^dx 
Jr 

where (x, y) T{x, y) is the integral kernej^ of a trace-class operator T considered as an 
integral operator on -L^(M); usually one denotes formally T{x,y) = {x\T\y). For future use, we 
define a linear space 

T^n) := {T G T{n) I f G S{M.^)}. 

Note that, the mapping T^{H) 3T ^ W'^ G 5(]R^) is bijective (see, e.g. prop. 4 of [2]). 
For any K and T G T{H) we get 

tr[TE^^(Z)] = - / ti[TD{z,)KD{z,y]d''z, = - [ pl{z,)d''z, 

JZ JZ 



where, by equations (2.4) and (2.1), the density 

pI[z^) := ii[TD{zi)KD{ziy]=A'Kfw^''^-''\z2)W'^''^'''>\z2)dh2 

Jc 



An I W''\zi/2 + Z2)W^{-Zi/2 + Z2)d'^Z2 = An / W^{z)W^{z - Zi)d'z 



c 



and Pxi^i) = Pt{~^i)- Note that p^ is a well-defined continuous function even if K is only a 



bounded operator. Thus, from equation (2.3), one has 



= ^p^u- 
Zn 

2.2. The POM associated to homodyne tomography. Let Q := [0, tt) x M so that gives 
a parametrization for the projective space (see. Appendix). Define the normalized positive 
operator measure (POM) E^t : B{Q) (^{'H) associated to homodyne tomography [2] as 

£^1(6 X i?) := - / Qe{R)dd = - [ e'^^Qo{R)e-'^^de, 6 G S([0,7r)), R G S(M). 

Obviously, one can define similarly a POM associated to the double covering x M of as 

Eht(e xR):=^ [ Q0{R)de = ^E^tiQ n [0,7r) x R) + lE^t{[{e n [7r,27r)) - n] x {-R)) 
Zn Jq Z Z 

where now 6 G -B([0, 2n)) and R G B{M). 



^That is, for all i/" G ^^(K), {T'4^){x) = J^f{x,y)^p{y)dy for dx-almost all a; G M. More precisely, T is an 
equivalence class of functions which satisfy this condition. 
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For any T G T(H), let p^^{9,r) be the density of the complex measure tr [TEht] with respect 
to {2-n)'^dedr [2j. For example, smce 

{n\E^,{Q X R)\m) = ^ [ e'^"^'de [ K{r)hm{r)dr, 

it follows that 

pt><"l(e,r) = e^("-)^/i„(r)/i^(r) = e^("-")^i><"l(0, r). 

The next proposition shows that, if T is regular enough, one can express p^^. using the matrix 
elements of T. 

Proposition 1. For any T G T^{T-[) we have: 

(1) the double sequence {Tmn)m,nm is rapidly decreasing^ 

(2) the integral kernel T{x,y) = J2m n=o'^-m.nhm{x)hn{y) where the double series converges 
absolutely for all x, y E ]R£| 

(3) plt{0,r) = YlZ,n=QTrnnhn{r)hm{r) = f{r,r) for (almost) all r G M, 

(4) the function r i— T(r, r) belongs to S{R), 

(5) e-^^^Te^^^ G T^{n) for all 6 eR, 

(6) pl{e,r)=pi;"'''^-'"'{0,r) = j:Zn=oTmne'^''-'^^'h^{r)hUr) for (almost) all 9 E [0,2n) 
and r G M, 

(7) the Wigner function W'^ G ^(M^); 

(8) the Radon transform {RW''^){6, r) = p'^^^iO, r) for (almost) all 9 G [0, 2tx) and r G M. 
Proof. (1) is a standard result, see e.g. theorem V.13]. Item (2) follows from the facts that 

su-p{\hn{x)\ >0|nGN, xgM}<1 

(see, equation (22.14.17) of [1] in p. 787) and X]mn=o \^mn\ < cxo; moreover, for all ip = 

c h Ic P < oo 



/ f{x,y)ilj{y)dy = ^ Tmnhm{x) / /i„(?/)V'(?/)(i?/ = ^ /^^(x) ^ T^^c^ = (TV')(a;) 

-^■^ m,n=0 '^'^ m=0 n=0 

for almost all x G M. Item (3) follows from 

tr [TQo{R)] = [ f{r,r)dr, R G B{R). 

^That is, the series X^m n=o "^^^"^'^1^™"!^ converge for all j, k eN. 

'^Obviously, here we have chosen a representative of the equivalence class of kernels T; if one chooses another 
representative then the results of this proposition holds only for almost all points. 



QUANTUM TOMOGRAPHY AND GENERALIZED MARKOV KERNELS 7 

Since, for any / G iS(]R^), the mapping r ^-^ f{r,r) belongs to iS(]R), item (4) follows. Item (5) 
is obvious and (6) is easy to see since Qe{R) = e'^^Qo{R)e-''^^ for all R G i3(M), 6 eR. Item 
(7) is a standard result (see, e.g. prop. 4 of [2]) and (8) is proved in p, prop. 2]. □ 

3. Results 

In this section, we will find different connections between phase space observables E^, K G 
T(H)^, and Eht, the observable associated with homodyne tomography. Our results are based 
on the properties of the Radon transform R and on the Hilbert transform H which are shortly 
introduced in the Appendix. The main mathematical results to be used are the convolution 



theorem (3.8), Radon inversion theorem (3.10), and Plancherel formula (3.11). 



3.1. Convolution theorem and Radon inversion theorem. Let T, K E '^{'H). The 

density of Z i— > tr [TExiZ)] at zi = (qi + ipi)/^/2 G C is 



p^(zi) := ti [T D (zi) K D (zi)*] = 4n W^{z)W^{z ~ zi)d^z 



rKi 

vv \z)yv 

IC 

= 27r / / W'^{q,p)W^{q - qi,p - pi)dqdp = 2-^ // {qi - q,pi - p)W'^{q,p)dqdp 
J Jr2 J Jr2 

= 27r{W^ * W^){qi,pi) = 27t{W^ * W^){qi,pi) 

where W^{q,p) := W^{—q, —p). If W'^ and are integrable then it follows from proposition 
2 of [2] that 

{RW'^){e,r)=pl{e,r) 

for (i^fir-almost all {6,r), and similarly for . Moreover, one gets from the convolution 



theorem (3.8) that 



(3.1) (Rp^)(^,r) =27r / pHO, s)p^,ie, s - r)ds 

for d^dr-almost all {9,r). Note that, formally, replacing]?^ with = ^Pm (which is only a 
quasiprobabihty distribution for non-Gaussian states) in the above equation, we get 

{RW'')(e,r)=pl{9,r)= [ pl{9, s)5{s - r)ds 

Jr 

and p™{6,r) = S{r) where S{s — r)ds is the Dirac measure concentrated on r, so that, the 
phase space distribution is a kind of a 'smoothed version' of quadature data (associated 
with W^). 
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Keeping equation (3.1) in mind, one gets from Radon inversion theorem (3.10) that 



(3.2) pI = ^R*A{Rpl) 



T K 

Pht' Pht 



if is a Scwartz function]^ Thus, the probabihty distribution where T, K E T^{7i) fl 
T{T-l)i tie constructed by using homodyne tomography data for both T and K. Physically, 
this means that to reconstruct the measurement outcome statistics pj^ of an 8-port homodyne 
detector, one needs to measure the statistics p^t separately by using balanced homodyne 



detection for quandratures, and then use formulas (3.1) and (3.2). In some cases, this could 
be easier in the sense that, if one needs to prepare or measure the kernel state K, one can use 
quadrature measurements (tomography) to get (and find the state K if it is unknown a 
priori). Then pj^ can be obtained directly either by using 8-port detection (since K is known) 
or by measuring p^^ via quantum tomography. Since that pj^ is (almost) symmetric in the 
interchange of K and T, that is, pj^{zi) = pff{—zi), mathematically, the signal state T and the 
kernel state K are in the same role in 8-port homodyne detection. 

3.2. Plancherel formula and Markov kernels. Let K G T{H) be such that G L^(]R^) 
and {q,p) G M^, (gi,Pi) G M^, and 9 G [0, 27r). If one denotes W^p^{q,p) := W^{q - qi,p-pi) 
and defines parameters a := qi cos 9 + pi sin 9, b := — gi sin 6* + pi cos 9, one gets 

i^^i^upi ) r) = [ W^^ p^ (r cos ^ - t sin ^, r sin ^ + t cos 9)dt 
Jr 

W^{r cos 9 — t sin 9 — qi, r sin 9 + t cos 9 — q2)dt 

/ W^{{r- a) cos 9 — {t — b) sin 9, (r — a) sin 9 + {t-b) cos9)d{t-b) 
Jr 

,r — a) = {RW^){9, r — qi cos 9 — pi sin 9) 
^^(6*, r — qi cos9 — pi sin 6'). 



Assume then that K,T e T^{n) and define M^f := A{RW^pJ (see, equation (g in the 
Appendix). Since 

= 2vr // W'^{q,p)W^{q-qi,p-pi)dqdp 



^This holds, e.g., when T, K e T^{n) since then is a convolution of Schwartz functions and 
and thus a Schwartz function (see, proposition jl| . 
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we get from Plancherel formula (3.11) that 

Jo JM 

Hence, for any Z e S(M^), 

tr[rEK(^)] = 7^ / pl{qi,Pi)dqidpi = ^ [ [ [ Mlp^{9,r)dtT[TEu{9,r)]dqidp, 
Jz Jz Jo Jr 

If Z is compact, it follows from Fubini's theorem, that one can change the order of the above 

integrations and get for operator measures Ex and Eht (by the density of T^iTi) C TiTi)) that 



2lT 



27r 



M^^^p^{e,r)dq^dpi 

z 



2lT 



dEUO,r)= / / MK{Z;e,r)dEu{e,r) 



--■.MK{Z;d,r) 



(weakly) for all compact Z C M^, where A4k is a generalized Markov kernel^ Since Aix is not 
necessarily a nonnegative function (see examples 1 and 2), AixiZ; 6,r) is not true conditional 
or transition probability. However, one could say that the measurement Ek is subordinate to 
the measurement Eht (see, [9]) or Ex is a postprocessing of Eht. 

Since pl^{qi,Pi) = Em,n.=o ^m"PiL>(n.|(^i'Pi) for all T G T(H) (where the double series 
converges absolutely) it follows that (weakly) 

/ MK{Z-e,r)dEu{e,r) = Ek(Z) = ir„„E|^)<„l(Z) 

Jo JK „_n 



m,n=0 



= ^ / j M\^)i^n\{Z]e,r)dE^^{e,r) 

m,n=0 ■^'^ 

for all compact Z C M^. Thus, one needs to find the density p_^{6,r) of the Markov kernel 
only in the case of i^' = \m){n\ where m, n G N. Before doing that we study some properties 
of the densities which are immediate from the definitions. 

We notice that, when gi = = pi, one gets P/^(0, 0) = tr [TK] where T, K E TiTi). Thus, 
if T, e T^iU), 

^ ^ ^ M^Je,r)pUe,r)dedr = tT[TK]. 



2vr JO 

Especially, by choosing K = \m){n\, m, G N, one recovers the matrix element T^m from the 
tomographic data p-[^ as follows: 

Tnm = :^ r [ Ml,-^^-\9,r)pl,{9,r)d9dr 

^TT Jo Ju 



%ote that the effects Eif (Z) of compact sets Z fully determine the POM Ek , K € T'^ n T(n)^ 



27r 



M^oi^,r)e''^^-''^^hi{r)hk{r)d9dr 
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(see, section 5 of [13] and references therein), li T = \k){l\, k, I & N, then 

Kik = ^ r [ M^^,{9,r)ptl^\9,r)d9dT = ^ 

since ^^^'^^\9^r) = e''^^~^^^hi{r)hk{r). Combining the above two results, we get the orthogonality 
condition: 

(3.3) ^ r [ AC^^''\9,r)e^^'-'^'hi{r)hk{r)d9dr = 5^,,5„,, 

^TT Jo 



(where 5m,i is the Kronecker delta). Next we calculate the structure of M^p_^{9,r) in the case 
of ii" = \m){n\. 

Recall that M^f = A{RW^^p^) where {RW^^pJ{9,r) = p^^{9, r-a) and a = cos9+pi sm9. 
Since P^-^'^^^^{9,r) = e^'^^~^^^hn{r)hm{r) one gets 



(RW^Si"')(^,^) = e^^"-")^/i„(r - a)hUr - a). 



Since the Hilbert transform H is shift invariant and commutes with derivation (see Appendix) 
one needs only to find H(/i„/im). 

Consider first the case n = m = (i.e. the Husimi Q-function case where K = |0)(0|). Then 



ho{x) 



--e and ho{xf 



1 

e . 



Let daw(r) := e~^^ e^^dx be the Dawson's integral which is related to the imaginary error 
function erfi in the obvious way (see, e.g. [U HH]). The result of the next lemma is well-known 
but, for completeness, we give an alternative proof for it. 



Lemma 1. For all r G 



[H{hl)]{r) = -e-"' / e'^'dx = -daw(? 



Proof. Let 



Since 



dfjr) 
dr 



f{r):= ml)]{r) 



lim 



7J-3/2 e^o+ _/ 2!6H r — X 

\r — x\> 



1 _ 2 

e erfi(r). 



dx 



lim 



—2x6 



I r — xl >e 



7r3/2 
2 



2^ f Q 



lim / dx — lim 

e-»0+ / xgR r — X e-^0+ J a:eE r — X 

r — x[>€ |r — a:|>e 



TT 



3/2 



e ^ dx — r lim 



e^0+ _/ i:6R r — X 

\r — x\'>e 



TT 



2r/(r) 
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and 



it follows that 



where c G M. But 



- — daw(r) = 2r— daw(r) 



2 

/(r) = — daw(r) + ce 

TT 



/(O) = li"^ / dx = — lim / dy = = + c 

|a;|>e 

SO that c = and the lemma follows. □ 

Since {RwS!fJ){9,r) = Mr-a)]^ (A^)(^,r) = 7r|:(H^,)(r), and = A{Rw!!f,), the 

above lemma implies that 

^liip? i0,r) = 2— daw(r - gi cos ^ - pi sin 6) = Y, 2H2k)\ ^^^^'' ~ """^ ^ ~ ^'"^ 

^ fc=0 ^ 

is an analytic function vanishing at infinity r ±oo (see, lemma 1 of pjj). Moreover, 

r.27r 



E|0>(0|(^) 



1 

2^ 



d^^Ae.r) 



(see a direct verification from ^I3j). Next we consider the general case when K = \m){n\ where 
m, n G N. 



Proposition 2. For all r G M, 



^^( ^- (-1)"^" Inlml"":^""^ 2" d"+""^^(2daw(r)) 

D=0 ^ ^ ^ ^ 

Proof. By the Feldheim's identity [TJ eq. (1.4)] 

min{n,m} ^t, 

Hn{r)Hm{r) = n\m\ — r:Hn+m-2vir 

' ' Diin — Dllim — 7)11 

and Rodrigue's formula 



v\[n — vjlim — vjl 



i/,(r) = (-l)V^'^" 



dr^ 

one can write 

1 _ 2 

hn{r)hm{r) = =H„{r)Hra{r)e 
V2"+'^n!m!7r 



mm{n,m} 

n!m! ^ 2" / ^ -r^ 



V2'^+™n!m!7r v\{n - v)\{m - v)\ 



1 / i rmin{n,m} jn+m-2«; -r^ 

1 /n!m! ^ 2^ . ^.„+^_2.^"+" '"e 
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Since e"''^ = ^/vr^o('")^ and H commutes with derivation, one gets 

n^nntim) [ i) y2„+^ ^ - v)\{m - v'y. rfr"+— 2- 

and the proposition follows from lemma [T] □ 

For any smooth function / : M ^ M we let f^^^ denote the sth derivative function of / (and 
/(o) := /). Let Y := 2daw^^\ As shown in the Appendix of [II], for all p G N, 

(3.4) r(^-)(r) = (-ir2-f:^^^^2.(r), 

(3.5) Yi^^-'\r) = (-l)^"^2-f:^^|^^^i^2..i(r) 

for all r G M. Finally, we are ready to calculate the density M]™pi"' of the kernel associated 
with \m){n\: from the previous proposition one gets 
Ml"'J<"l(e,r) = ^c'^n~m)e 9[^{hnhni)]{r - a) 



I j j-min{n,m} 

^ ' ^ ' V 2^^+"^ ^ v\{n - v)\im - v)\ ^ ' 

where a = q\ cos ^ + pi sin Q. Immediately one sees that M 9 r M^^i^'^\^ , r) G M is analytic 
and vanishes at infinity r ±00. Moreover, 



^Si"' (^, r) = e'("— )ejvf l-)<"l (0, r - gi cos e - pi sin 9) = M'^p^^' (^, r) 



and the orthogonality relation (3.3) reduces to 



1 f 

Mff ^"'(0,r)/iKr)/^fc(r)rfr — / e^C— = 5„,i5n,k. 



^ ^ 



— Sl-m-k + n,0 

Obviously, without restricting generality, we may assume that 6' = 0, gi = 0, and pi = in the 
next theorem: 

Theorem 1. If n + m is even, then 

^^^0,0 lU,rj- / , , / , okfou\} n TT — , vTi ^2k[r), 

k=it^m)/2^^^^y- {k--,{n + m))\ 

and if n + m is odd, then 

M\-)i^\(nr) (-1)^"^"^"'^^' (-1)^ {k + Un-m + l))\{k + l{m-n + l))\ 

^^^0,0 lU,rj- , — > ofc/oL I i\| 71 TT — , TvTj ^2k+i[r) 

^^^^^_^^^^2'=(2A: + 1)! (/, _ + ^ _ 1)) 1 

for all r G M. 
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Proof. Suppose that n + m is even so that one can write n + m = 2/ where / G N. Then from 



equations (3.6) and (3.4) one gets 



" min{n,m} 



y 2n+m yllji — Vjlim — V)\ 



min{n,m} 



\/n\m\ 

2' ^ t;!(n-f)!(m-t;)!^ " ^ 2'=(2fc)! 



;^2'=(2A;)! t;!(n-t;)!(m-f)! 

k=0 ^ ' v=0 \ / \ 



n\ fk + l — v 

m — V 



(see eq. (5) of [18] in p. 

1^^V^ (-1)' f7 , ,(A; + /-m)!(A; + /-n)! 

Vn m (-1) > , . , ,. H2kir)—, — ry, 

V ^Z^2'=(2A;)! n!m!(fc + / - n - m)! 



k=l 



-ly ^ (-1)^ (A; + /-m)!(A; + /-n)! 



^/^^2'=(2A;)! (A; + /-n-m)! 

and the first equation follows by substituting / = (n + m)/2 into the above equation. 
Similarly, if n + m = 2/ + 1, / G N, then from (3.6) and ( |3.5 ) it follows that 



ain{n,m} 



~\ T ^ ( T\k min{n,m} , , , , -.x, 

n!m! ^ (—1) rr [k + 1 — v + ly. 

"2~^2'=(2A; + 1)! ^ t;!(n - t;)!(m - t;)! 

iN^Y^ (-1)^ rr , ,(A; + l + /-m)! + l + 

( m ) 

i-iy ^ (-1)^ ik + l + l-m)\{k + l + l-n)\ 
" v/27^^2*^(2fc + l)! (A; + l + Z-n-m)! 

and the second equation has been proved. □ 



Remark 1. Since Y is an analytic function (see, eq. (A. 6) of [TT]) one sees from equation (3.6) 
that r^M'^ <"l(0 , r) is an analytic function. Indeed, its Maclaurin series is easy to compute 
by using the series 

(3-7) nr)=2j:^-^i2rr 
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which imphes that, for all p G N, 



2u^2u~p 



u>p/2 



For example, if n + m = 2/, / G N, one gets 



" min{n,m} 



" min{n,m} 



22' ^ (2M-2/ + 2t;)! 

i=0 ^ ^ i)=0 ^ ^ ^ ^ 

3.3. Examples. 

Example 1. Let us consider the case when K = where is a Schrodinger cat state: 

^ = ;^(|0) + |i))- 



From Q we get Ml°i^°l(^,r) = r(r), Ml°i^^l(^,r) = -e^^rW(r)/V2, and Ml^i<^l(e,r) 



Y{r) + Y^'^\r)/2, so that 

Ml5<'^l(^,r) = Y{r)- ^Y^^\r) cose +-Y^^\r) 

V 2 4 

= 2(r + v/2cos^)[r + (l-2r2)daw(r)]. 

The function [0, vr) x M 9 {6, r) i— * Mq'^^^'^' (^, r) G M is plotted in the next picture by interpreting 
{9,r) as the polar coordinates of M?. 

Example 2. In this example, we study the Cahill-Glauber s-parametrized quasiprobability 
distributions [3]. Instead of using s as a parameter, we define A := (s + l)/(s — 1) and assume 
that A G [-1,1) (that is, s = (A + 1)/(A-1) G (-oo,0]). Let T G T{n). Now the s-distribution 
isp]^^ where the kernel Kx := (1-A) E^=o^"l^)(^l l^, eq. (6.23)]. Note that Kx G T^iH) only 
when A G (—1, 1), and then tri^A = 1- For A = (and s = —1), Kq = |0)(0|, and pj^^ is the 
Husimi Q-function. When A = — 1 (and s = 0) we have K^i = 211 and = P2u = 2ttW'^ . 
Note that, for A G (—1, 1), one can define an operator measure E;^^ and calculate MQ^iO^r) as 
before, but when A = — 1 the corresponding E/^_^ is only a generalized operator measure, that 
is, a sesquilinear form valued measure [H]. Hence, we assume that A G (—1, 1) and finally take 
the limit A ^ — 1. 
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(a) 



(b) 



Figure 1. The density mI^q^'^^O, r) of the generahzed Markov kernel in the case 
of a Schrodinger cat state ip = '^{\^) + l^))- 



By theorem [TJ one sees that 



r) = (l-A)5^A"Ml:><"l(^,r) 



n=0 

oo 



— fc — 71, 



k=0 



n=0 



n\{k — n)\ 



k=0 



It is easy to calculate the Maclaurin series of r i-^ M^q{6, r) (see eq. (A.l) and the text above 
eq. (A. 6) of [II]): we get 

u+l 



-1/. 



From equation (3.7) one sees that 



Mo^o^(^,r) = —Y{r/V^s) = — M^^l^'^ {6 , r / V^s) . 



When A ^ -1, s ^ 0, 



so that the density of the Markov kernel is not defined in the Wigner function limit, as expected. 
Note that does not define the Dirac (5-distribution in the limit (even though it is an 
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increasing peak, see the next picture), since J^Y{r)dr = implies that 

Jr 

for all A G (—1, 1). However, for any G one gets 

but cannot change the order of taking the limit and integration. 



Figure 2. The density r i— > MQQ{9,r) of the generalized Markov kernel (for all 
9 G [0, 27r)). The higher curve corresponds to the case s = and for the lower 
curve s = —1. 

Example 3. In this example, we show that the generalized Markov kernel does not exist for all 
states. Choose Kr/ = \Ti){ri\ where t] = Xlo,i] ^ -Z^^(M) (the characteristic function of [0, 1]) is a 
unit vector. The first observations are that Krj ^ T'^{T-C) and the Wigner function W^"^ of Kj^ is 
not integrable. Thus, the equation {RW^^){9,r) = p^l\9,r) is not necessarily valid. However, 
we can calculate p^^{9,r). For example, p^''(0,r) = [x[o,i](^)]^ = X[o,i](^)- Although the Radon 
inversion theorem and Plancherel formula obviously fail, we can try to derive M^^^(6',r) as 
before we did for T'^(7i)-kernels. But now dx[o,i\{x) / dx = 6{x) — 6{x — 1) must be understood 
as a (tempered) distribution and a formal calculation shows that 

which is singular at r = and r = 1. On the other hand, 

7r{Hx[o,i]){r) = In 

and a we get the same result: 
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It is interesting to note that E^^ is not informationally complete, since the support of 

{q,p)^tT[K,D{q,p)] = {r]\D{q,p)r]) = 6-'"^/^ f e'^^^Xm^x - q)dx 

Jo 

is not the whole [Sj [12] . It is an open question how the informationally completeness of a 
phase space observable is connected with the existence of the generalized Markov kernel 
Mk- 

Appendix: Radon and Hilbert transforms 

The basic properties of the Radon and Hilbert transforms, R and H, respectively, are collected 
in this Appendix. There is plenty of literature on these transforms. For example, [8] and 
[6] essentially contain the results for R whereas and [15] are good references for H. The 
elementary results of this Appendix which are easy to prove (e.g., by changing integration 
variables and by using Fubini's theorem), are denoted by {-k). 

Let be the set of straight lines in M^. Any line C,{0,r) := {{x,y) G | xcos6 + ysin6 = r} 
can be parametrized with two real numbers 6 G [0,27r) and r G M, and becomes a smooth 
two-dimensional manifold. Since ^{—6,—r) = ^{6,r), P^ has a double covering x M ^ P^ 
where is the unit circle {(cos 6*, sin 0) G | G [0,27r)}. An unambiguous parametrization 
of P^ is then given with the same parameters as above but restricting their domains to be 
either 6 G [0, vr) and r G M, or G [0, 27r) and r > 0. From now on we use the double 
covering space §^ x M instead of P^ and, for example, identify functions P^ ^ C with functions 
/ : [0,27r) X M ^ C with the property /(-0, -r) = f{e,r) for all 6 G [0, 27r) and r G M. 
Moreover, without any further mention, we equip §^ x M with the measure {27i)~^d6 dr. 

Let iS(P^) rapidly decreasing smooth functions P^ ^ C (obviously in the direction of growing 
r). Let iSiy(P^) consist of functions G 5(P^) such that, for all k G Z+, the mapping 6 \—>- 
J^(f{6,r)r''dr is a homogeneous polynomial of degree k with respect to 'variables' cos^ and 
sin 6'. 

The Radon transformation is a continuous (^) linear mapping R : L^(]R^) L^{S^ x M) given 

by 

(R/)(6',r) := / /(rcos6' - tsin6', r sin6' + tcos6')rft 
Jr 

for dOdr-almost all {6, r); here we consider the above L^-spaces as Banach spaces equipped with 
the stardard L^-norms. We have the following Schwartz theorem [HI theorem 2.4]: 

R|5(M2) : 5(M^) Sh{P^) is bijective. 
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The operator R has an adjoint operator, the so-called backprojection operator, R* defined as 

{R*(f){x,y) := I (p{6 , X cos 9 + y sin 6) — 
Jo 27r 

where, e.g., ip G L^{S^ x M). Indeed, for suitable classes of functions / : ^ C and 

{p : X R ^ C one can easily verify (^) that (R/|(y9)j;^2(§ixK) = {f\R*ip)L2(^ 

Let f,g^ L-^{M.'^), and define the convolution f * g G L-^{M.'^) as 



if * 9){x,y) := // f{x -t,y -u)g{t,u)dtdu. 

J JM? 

We have the following convolution theorem (*): 

(3.8) [R(/*^)](^,r)= [{Rf){e,r-s){Rg){e,s)ds. 

Define the Hilbert transform H which is a bounded linear operator on L^(]R) given by 
(H^.)(r) := i lim / = -i 1„„ f ijt^vl^^t^,, 

TT e^0+ J xsR r — X TT e-'0+ J ^ y 

I r — X I > e 

for almost all r G M and the above limit exists also in the L -norm H theorem 8.1.7]|3 The 
Hilbert transform has the following properties: 

(1) H is shift invariant that is, 

(HV'a)(r) = (H^)(r-a), a, r G M, 

where if) G L^(M) and ipa{x) '■= ip{x — a) for all a; G M. 

(2) H commutes with derivation [U prop. 8.3.7], that is, 

l_l . d'ip\ dHxj) 



dx J dr ^ 

where ip G L^(]R) is absolutely continuous and dip/dx G L^(]R) (that is, when belong 
to the domain of the momentum operator P). 

For any Lp G 5h(P^) define (pe{x) := ip{9,x) and 

(3.9) (A^)(0,r):=.H(!^)(r)= hm / ^^?f^dx = ^(H^.Xr). 

\ dx J e^o+ J xm r — X ox dr 



r — a: I >e 



p2N 



We have the Radon inversion theorem P, theorem 3.6]: for all / G 5( 
(3.10) / = ;^R*AR/, 

Note that some authors multiply H by the imaginary unit i in the definition of the Hilbert transform to 
get a self-adjoint operator. 



QUANTUM TOMOGRAPHY AND GENERALIZED MARKOV KERNELS 19 
which imphes the following Plancherel formula: 

II f{x,y)g{x,y)dxdy = (/|^)i2(]R2) = ;^(/|R*AR^)i2(K2) = -^(R/|AR^)i2(gixK) 

J J^2 ZTT Zn 

(3.11) = j^ j^^jRfWy){^^9){0^r)dddr 

for all f,ge ^(M^). 
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